Gauge Field and Confinement-Deconfinement Transition in Hydrogen-bonded ferroelectrics 
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Quantum melting of ferroelectric moment in the frustrated hydrogen-bonded system with "ice rule" is studied 
theoretically by using the quantum Monte Carlo simulation. The large number of nearly degenerate configura- 
tions are described as the gauge degrees of freedom, i.e., the model is mapped to a lattice gauge theory which 
shows the confinement-deconfinment transition (CDT). The dipole-dipole interaction J2, on the other hand, 
explicitly breaks the gauge symmetry leading to the ferroelectric transition (FT). It is found that the crossover 
from FT to CDT manifests itself in the reduced correlation length of the polarization £ft ~ A(K — K c )~ v 
with A oc \/~J% while K c and v remains finite in the limit J2 — > 0. In contrast, the susceptibility \ ar, d the 
spontaneous polarization continues smoothly in this limit. Moreover, CDT is related to the shape or volume of 
the molecules, and the its correlation length £cdt is not reduced in the limit of J2 — > 0. A proposal for the 
neutron scattering or X-ray and Raman experiments is given to confirm this prediction. 

PACS numbers: 64.70.K-, 64.60.-i,11.15.-q 



The hydrogen-bonded systems are one of the ideal laborato- 
ries to study the quantum tunneling. Especially, the ferroelec- 
tric properties of these systems attract much attention since the 
old work by Slater on KDP HI. The quantum melting of the 
ferroelectric order to result in the quantum paraelectricity is 
a rather common phenomenon observed in several hydrogen- 
bonded ferroelectrics JIJQ, which is usually described by the 
transverse Ising model as given by 

H=-Y d Jii<rtv z j -KY,° x , (1) 

ij i 

where a z = ±1 specifies the two positions of the hydrogen 
atoms, Jij is the dipole-dipole interaction, and K represents 
the tunnelling matrix element. These two interactions com- 
pete with each other, and by increasing K, there occurs a 
phase transition from the ordered state to the quantum dis- 
ordered phase. It is well known that the model Eq.([T]i in 
rf-spatial dimensions at zero temperature is equivalent to the 
(d+ l)-dimensional classical Ising model at finite temperature 
T with K playing the role of T EE). 

On the other hand, it often happens that the constraints 
are significant to the hydrogen-bonded systems. Actually, the 
hydrogen positions in the representative system KDP are al- 
ready subject to the constraint, i.e., so called "ice rule" (TJ. 
Namely, only two of the four hydrogen atoms next to a tetra- 
hedron are approaching to the center for the low energy sec- 
tor. Similar constraint is also relevant to the recently studied 
quasi-two dimensional antiferroelectric squaric acid, where 
the square molecule is surrounded by 4 molecules with hy- 
drogen bonds J5], and the two-in-two-out configurations are 
energetically stable. This "ice rule" is the generalization of 
the hydrogen bonds in ice leading to the macroscopic degener- 
acy of the ground state configurations as discussed by Pauling 
long time ago (8). Therefore, a keen issue is how this macro- 
scopic degeneracy of the low energy states affects the physical 
properties of the hydrogen bonded systems, which we address 
in this paper. 



The constraints imposed on the physical variables are more 
common phenomenon found in many other cases. Frustrated 
magnets are one of such examples, where some of the macro- 
scopically degenerate spin configurations are selected as the 
lowest energy states. Spin ice in pyrochlore ferromagnet is 
a representative example, in which the hydrogen position is 
replaced by the direction of the spin, and the "ice rule" is ap- 
plied here also. This fact leads to an interesting phenomena, 
e.g., the absence of the long range ordering down to zero tem- 
perature and the deconfined magnetic monopoles as the ex- 
citations J5]- These are described well by the gauge theory 
representing the constraints within the framework of the clas- 
sical statistical mechanics. Quantum effects on the spin ice 
model have been attracting intense interests recently ifTUl [TO . 

In this paper, we study theoretically the phase transition 
in a (2+l)-dimensional model of the hydrogen-bonded ferro- 
electrics with the ice rule. This model shows two types of 
phase transitions, i.e., the confinement-deconfinement transi- 
tion (CDT) and the ferroelectric transition (FT), and the issue 
is how these two are connected. As discussed below, the latter 
is induced by the dipole-dipole interaction J 2 which explicitly 
breaks the gauge symmetry, while the former exists with the 
pure gauge model. 

Inspired by the squaric acid, we consider the following 
model 

H = H Q + H 1 +H 2 , (2) 

where 

H a =-J ^afai < 4ai-K^2af, (3) 

□ t 

Hi = Ji^K^+a 2 X), (4) 

□ 

H 2 =-J 2 Pa-Pb- (5) 

<AB> 

in the lattice given in Fig. [TJ\. The summation of □ in Eq. ([3]) 
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and |4]) is over the blue plaquetts which resemble the H2SQ 
molecules 0. The sites are labeled by numbers and the con- 
vention is given in Fig. [TJ\. On each site, there is a hydrogen 
shared by two neighboring molecules, representing the hydro- 
gen bond. a z is used to denote the motion of protons: If a 
hydrogen is closer to the plaquette of group A, it is defined 
as the "+" state, otherwise it is a "— " state. The summation 
in Eq. ^ is over the nearest-neighbor plaquettes, and P, are 
defined as the following 
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Pa* = 4K 



<7n — <7o — 



Pa v = 7 (erf + °3 - °i - a %) 



for the plaquettes of group A, and 
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Pbx = 7(03 + 0-4 



— CTi — a 



2)- 



Pbij = 7 (of + 0% - cr| - 0-3) 



(6) 
(7) 

(8) 
(9) 



for the plaguettes of group B, representing the polarization 
vectors. 

The model described by H only is nothing but the quan- 
tum Hamiltonian of the (2+l)-dimensional Ising gauge theory, 
where a z is the bond variable of the dual lattice illustrated by 
the dash plaguette in Fig. [T}\. It is well known that it is dual 
to a three-dimensional classical Ising model. When K = 0, 
it becomes a two-dimensional Ising gauge theory dual to the 
one-dimensional Ising model J6] |7|. There is macroscopic 
ground state degeneracy. In each plaquette, there are 8 dif- 
ferent configurations in the low energy sector illustrated in 
Fig. [TJ3. The addition of Hi lifts the degeneracy so that the 
states of (e) to (h) are favored. States (e) to (h) are particularly 
interesting because they carry finite dipole moments. The di- 
rection of the dipole moments for the plaquettes of group A 
are shown in red arrows in Fig . [Tj3 . Finally, a nearest-neighbor 
dipole interaction is introduced by H2. 

Let us consider the classical model first, namely K = 0. 
The finite-temperature phase diagram is summarized in Fig. [2] 
For Jx = J 2 = 0, the system is always disordered, since the 
model is equivalent to a one-dimensional Ising model by an 
appropriate gauge choice, i.e., there is no phase transition at 
finite temperature. For J% >> Jo and J 2 = 0, the gauge sym- 
metry is partially broken. Namely, the model becomes infi- 
nite number of decoupled one-dimensional antiferromagnetic 
Ising chains, and is again equivalent to the 2D Ising gauge 
model. Therefore, the system is still disordered without finite- 
temperature phase transition. This can be explicitly shown by 
introducing the 77 variables defined in the dual lattice, for ex- 
ample erf = r)iT)j as shown in Fig.[TJ\. Then, to minimize the 
energy for Hq, the Hamiltonian H at J 2 = is given by 



H=-J n 



ViVimiih 



(10) 



where □ is the spatial plaquettes in the dual lattice. The effec- 
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tive Hamiltonian of Eq. ( 10 1 is a classical 2-dimensional Ising 



FIG. 1: (A) The squaric acid is layered structure, and H2SQ 
molecules span a two-dimensional network marked in blue. The 
black dots denote the hydrogen sites. We divide them into two groups 
to define the dipole moment. If the hydrogen is closer to (away from) 
the molecules of group A (B), it is a + state, otherwise it is a — state. 
The dual lattice is drawn in dash lines. (B) The ground state configu- 
ration of each plaquette in the model of Ho only. There are no dipole 
moment in the states (a) to (d), and there are finite dipole moments in 
the states (e) to (h). The directions of the dipole moment are drawn 
in red for the molecules of group A. 



gauge theory as mentioned above. Therefore, for J 2 = and 
K = 0, it is in the disordered phase over all temperature range 
shown as the pink region in Fig. [2] For J2 7^ 0, H2 introduces 
a two-dimensional Ising model, which breaks explicitly the 
gauge invariance. There is a ferroelectric phase transition at 
finite temperature. For small J2, the transition temperature T c 
behaves like ~ Ji/log( J\/ J 2 ) for Ji >> J , which vanishes 
as J2 goes to zero. 

Now let us turn to the quantum model. For K 7^ 0, the 
zero-temperature phase diagram can be summarized in Fig. [3] 
When Ji — J2 = 0, there is a confinement-deconfinement 
transition (CDT) at critical K c . Since the model is dual to 
a (2+l)-dimensional Ising model, the CDT is a second or- 
der phase transition SO- In order to demonstrate whether 
or not the CDT extends to finite J\ region (at J2 = 0), we 
perform the quantum Monte Carlo calculation in addition to 
the analysis using the r\ variables. The numerical results, or- 
ganised in the Supplementary Information, indicate that the 
second-order quantum phase transition extends to finite J\ re- 
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Ji/(Jo+ J,) 



Jl/(Jo+J2) 

FIG. 2: Finite temperature phase diagram of the classical model for 
K = 0. The disorder phase is drawn in pink, and the ferroelectric 
phase is drawn in green. Ji is the intra-molecular coupling, and J2 
is the dipole-dipole interaction. 



gion, which is consistent with the ^-variable analysis. The 
consistency is not surprising since the model behaves like a 
one-dimensional transverse-field Ising model in the J 1 — >• 00 
limit, where there is a second-order quantum phase transition 
separating the phase with critical correlation and the classical 
paramagnetic phase. The CDT extends to the finite J\ region 
and continues smoothly to the quantum phase transition at 
J 1 =00. As a short summary shown in Fig. [3] on the J 2 = 
plane, the phase diagram is divided into two regions. One is 
the deconfined phase and the other is the confined phase. 





deconfinement 
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FIG. 3: The quantum phase diagram of Eq. |2). For J2 = 0, there are 
two phases separated by a confinement-deconfinement phase transi- 
tion. For finite J2, the phase space is divided by a second-order 
ferroelectric phase transition. 

Introducing the dipole-dipole interaction J 2 7^ 0, the 
ground state develops a spontaneous polarization at small 
transverse field K. A quantum phase transition to the para- 



electric state can be driven by K. In Fig. [4] the dielec- 
tric constant is computed for 5 different values of J 2 . In 
Fig. [5] we show the results of the correlation length £ft for 
the dipole moments for K > K c . The ferroelectric transi- 
tion is a second-order phase transition because both dielec- 
tric constant and the correlation length diverge at K = K c . 
Therefore, the phase space is divided into two regions for fi- 
nite J 2 : a ferroelectric phase and a paraelectric phase. Ex- 
tending to the finite J 2 region, the deconfined phase at the 
zero- J 2 plane becomes the ferroelectric phase, and the con- 
fined phase becomes the paraelectric phase as shown in Fig. [3] 
Due to these non-trivial connections, how does the criticality 
of the confinement-deconfinement phase transition (CDT) of 
the gauge field affect the criticality of the ferroelectric phase 
transition (FT) is the main scope of this paper. 
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FIG. 4: The inverse of the dielectric susceptibility x f° r Ji = 
0.2,0.4,0.8, 1.2,2.0 is computed for J = 1 and Ji = 0.2. The 
dielectric susceptibility satisfies well the Curie-Weiss's law x = 
C I (K — K c ). In the inset, we extract C and K c as function of Ji. 
The curve in panel C uses the left y— axis, and the one for K c uses 
the right y— axis. C is almost 1 independent of J2, indicating the 
robustness of the ferroelectric transition for any finite J2. K c has a 
linear relation with J2, terminating at K c = 0.67 for J2 = 0, which 
is the K c for the confinement-deconfinement phase transition. 

As shown in Fig. [4] the dielectric susceptibility satisfies the 
Curie-Weiss's law X = C(K - K c )~ x with C = 1 for all 
J2, where C carries the information of the square of the po- 
larization moment. Our results indicate that the ferroelectric 
transition is robust for any finite J 2 , i.e, K c converges to a 
finite value. As shown in Fig. |5]\, the correlation length £ 
obeys £ FT — ( K ^ K y ver y nicely with v = 0.46, 0.42, 0.41, 
0.46, 0.55 for J 2 = 0.2, 0.4, 0.8, 1.2, 2.0, remaining finite as 
approaching to J 2 = 0. We believe that the fluctuation of the 
v comes from error bars in the estimation. Using the field- 
theoretical technique, detailed in the Supplementary Informa- 
tion, we obtain 



A 



J-2 



2J 2 + 4(2J - Jx)+X 



(11) 
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where A is the Lagrange multiplier to impose the constraint of 
the fixed number of spin states. Our numerical results can be 
fitted well with the analytical result by a proportional constant 
2.2 as shown in Fig . |5j3 . As the correlation length defines the 
critical region, our results indicate that the gauge field regu- 
lates the critical region as the system approaches to the critical 
point. Although the ferroelectric transition is robust for any fi- 
nite J2, the correlation vanishes as the system approaches to 



the CDT, consistent with Eq. ( 10 1 that any quantity which is 
not gauge invariant has the zero expectation value [6, 7|. This 
profound property provides a good measure of distance for a 
ferroelectric system in the vicinity of the CDT. The measure- 
ment of the correlation length by neutron scattering or X-ray 
scattering lfl2l \l3l toward the phase transition point can pro- 
vide a strong evidence that the system is near the CDT, i.e., 
the degree of frustration, in comparison with the susceptibil- 
ity and spontaneous polarization, both of which do not show 
any special behaviour. 
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FIG. 5: (A) The log-log plot of the correlation length £. v are ob- 
tained as 0.55, 0.46, 0.41, 0.42, 0.46 for J 2 = 2.0, 1.2, 0.8, 0.4, and 
0.2 respectively. As the dipole-dipole interaction J2 is switched on, 
the system immediately becomes an ordinary ferroelectric system. 
(B) The coefficient A of the correlation length as the function of J2. 
The blue-dot line is the fitting of the analytical result. 

Next let us turn to the CDT. The order parameter for CDT is 



intricate. It is a conventional local order parameter in the dual 
variable but a nonlocal order parameter in original cr's. For 
simplicity, we consider below the case of J\ = 0. The product 
of four cr z 's around a plaquette i, i.e., p t = a'^a^a^a^ is 
the gauge invariant and relevant quantity for CDT. Usually, 
the Wilson loop W(C) — fliPi' where the product over i 
runs over the all plaquettes inside the loop C, is the "order 
parameter" for CDT |6l|7). However, this quantity is rather 
difficult to be related to a physical observable in the present 
case. Instead, we consider the correlation of p, 's which shows 
the scaling behavior in the critical region as 



< pipj >~ R 



-(d-a/v) 



5(-R/£cdt), 



(12) 



where R = \Ri — Rj\ is the distance between the two pla- 
quettes i and j, while a and v are the critical exponent for 
the specific heat and the correlation length, respectively. The 
function g is a scaling function and £cdt is the correlation 
length of the dual 3D Ising model, which diverges as K ap- 
proaches to K c [7 |. We expect the same critical exponent v 
for £cdt as £ft for the FT. However, £cdt is not reduced in the 
limit of J2 — > 0, which means that we have two correlation 
lengths £ft and £cdt for small J2 . 

Now the question is how to detect £cdt experimentally. It 
is noted here that pi is concerned about the shape or volume 
of the molecule i. It is known that the change in the molecular 
shape results in the change in the frequencies of some molec- 
ular vibrations [12]- For simplicity, we consider the model 
where frequency of a molecular vibration flj depends on the 
plaquette variable as Hi = Slo + api. Then, we consider the 
Raman scattering spectrum due to this vibration mode, whose 
intensity I(q, lu) is given by 



i,3 



~iq-(R.i-R :j ) e i(Lu(t-t')-(l i t+n j t') > 



(13) 



where q is the momentum transfer between the incident and 
emitted lights and u> is the Stokes shift. Assuming the Gaus- 
sian fluctuation for p/s in the coarse grained model, we can 
obtain 



-^>exp[- 2 , 
L a 2 (l+<p iPj >) 

(14) 

Therefore, the width T q of the peak around oj = Oq depends 
on the momentum q as T q = r(q£cDT), which increases about 
a factor of \/2 when q increases across q ~ £cdt- ^ n ^ s 1 
dependence of the Raman spectrum is expected to give the 
information on £cdt, which should be different from £ft when 
the system is near the gauge symmetric model, i.e., the limit 
of large degree of degeneracy. 

In conclusion, the effect of the gauge field in the hydrogen- 
bonded ferroelectrics is very subtle. In the vicinity of the 
CDT, the gauge field suppresses the growth of the critical re- 
gion by regulating the correlation length. Our theory provides 
a scheme to uncover the shadow of the gauge field as well 
as to realise the accompanying CDT by identifying the two 
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length scales £ft and £cdt near the ferroelectric phase transi- 
tion. Finally, our theory may shed light on a possible gauge 
structure in the real ice, ice XI, the ferroelectric phase of ice. 
It is not only important in laboratories but also with astrophys- 
ical implications lfl4l[T5l . A future research direction toward 
the understanding of the gauge structure in ice is needed to be 
explored. 
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